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The exploration of the strong-interaction matter under extreme conditions is one of 
the main goals of relativistic heavy-ion collisions. We provide some of the main results 
on the novel properties of quark-gluon plasma, with particular focus given to the strong 
collectivity and the color opaqueness exhibited by such hot and dense matter produced 
in high-energy nuclear collisions at RHIC and the LHC. 
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1. Introduction 

The study of QCD matter under extreme conditions of temperatures and densi¬ 
ties is of great interest in the field of high-energy nuclear physics. The simulations 
from lattice QCD have demonstrated that as the temperature of the QCD matter 
is increased above a certain value Tc ~ 155 MeV™^ quarks and gluons that are 
confined in normal hadronic matter will be liberated, and form a novel state of 
matter which is usually referred to as quark-gluon plasma (QGP). QGP is believed 
to have existed at a few micro-microseconds after the Big Bang, when the temper¬ 
ature of the Universe was so high that quarks and gluons could not bind together. 
Such hot and dense nuclear matter may be created in the laboratories by colliding 
two heavy nuclei at ultra-relativistic energies, such as those performed at the Rela¬ 
tivistic Heavy-Ion Collider (RHIC) at Bookhaven National Laboratory (BNL) and 
the Large Hadron Collider (LHC) at European Organization for Nuclear Research 
(CERN). The main goal of high-energy heavy-ion experiments is to explore various 
novel properties of the hot and dense QGP, and map out the phase structure of the 
strong-interaction matter. With the center-of-mass collision energies up to 200 GeV 
per nucleon pair at RHIC and more than one order of magnitude higher at the LHC, 
the highest temperatures of the hot matter created in these experiments can reach 
^ 370 — 470 MeV,®^ which is well above the transition temperature predicted by 
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Fig. 1. (Color online) Illustration of different dynamical evolution stages in a typical heavy-ion 
collision. 


the lattice QCD calculationsPI^ 

Our current understanding of a typical heavy ion collision is schematically illus¬ 
trated in Fig. One starts from two fast-moving nuclei heading on each other. In 
the laboratory frame, the two colliding nuclei look like pancakes due to the Lorentz 
contraction. After two nuclei penetrate through each other, a large amount of en¬ 
ergy is deposited in the central collision zone and a highly-excited nuclear matter is 
produced. The produced matter experiences a pre-equilibrium evolution stage, and 
then reaches (a certain degree of) local thermal equilibrium. The equilibrated hot 
and dense QGP matter then expands like a relativistic fluid and cools down. When 
the temperature/density of the fluid drops below a certain value, it undergoes a 
transition from QGP to hadronic phase. The produced hadrons continue to inter¬ 
act with each other until the kinetic freezout, and then fly to the detectors. One can 
see that the deconfined color degrees of freedom in the QGP matter produced in 
energetic heavy-ion collisions are not directly observable, but manifest in the final 
state hadrons and other colorless objects. The essential task in the study of hot and 
dense QGP using relativistic heavy-ion collisions is to find clear and unambiguous 
connections between the produced QGP and the final (hadronic) observables, and 
to identify reliable signatures for the formation of QGP. 

Experimental results have shown that the QGP matter produced in these en¬ 
ergetic nuclear collisions exhibits many remarkable properties that differ from any¬ 
thing that has been seen before.l^lli^ On one hand, such hot and dense nuclear 
matter shows strong (anisotropic) collective behavior. The successfulness of rela¬ 
tivistic viscous hydrodynamics in the description of the space-time evolution of the 
bulk matter, especially the observed anisotropic collective flow, indicates that the 
produced QGP matter is strongly interacting and has very small specific shear vis¬ 
cosity. On the other hand, the produced matter shows strong opaqueness to the 
attenuation of high energy partons; this is usually called jet quenching. Model cal¬ 
culations based on parton energy loss which originates from the interaction between 
the propagating hard parton and the hot and dense QCD medium can explain quite 
well many phenomena associated with hard jets. These and other exciting findings 
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from RHIC and the LHC indicate that the studies of strong-interaction matter in 
high energy nuclear collisions have entered a completely new era. 

In this report, we mainly focus on two most important phenomena as men¬ 
tioned above, i.e., anisotropic collective flow and jet quenching in relativistic nu¬ 
clear collisions. In particular, we will discuss the recent efforts on the quantitative 
determination of two important transport coefficients, namely, the shear viscosity 
to entropy ratio (the specific shear viscosity) yy/s and the scaled jet quenching pa¬ 
rameter g/T^ (or q/s), via systematic phenomenological studies of soft particles 
produced from the bulk matter and high transverse momentum particles originat¬ 
ing from hard partonic jets at RHIC and the LHC. The precise determination of 
these two transport coefficients and their temperature dependences can provide us 
a quantitative estimation of the interaction nature of the hot and dense nuclear 
matter produced in heavy-ion collisions, e.g., how a weakly-coupled quark-gluon 
matter at sufficiently high temperatures turns into a strongly-coupled relativistic 
fluid at the energies achieved at RHIC and the LHC experiments.!^ 


2. Collectivity and Fluidity of QGP 

In a typical non-central nucleus-nucleus (A-A) collision, the produced hot and dense 
nuclear matter is anisotropic in the plane transverse to the bean direction (see Fig. 
[^. This geometric anisotropy is usually quantified by the second-order eccentricity 
£2 = — 2 ;^)/(j/^ -I- x^). Due to the interaction among the medium constituents, 

the initial geometric anisotropy will be converted into the anisotropy of the final 
state particle momentum distribution, which is usually quantified by the elliptic 
flow coefficient V 2 = {pi—p^)/{pi -l-Py).!^^^^The observed large elliptic anisotropy 
and the dynamical evolution of the hot and dense QGP matter have been well 
described by the simulations from relativistic hydrodynamics.!22HlZl This has led to 
the conclusion that at the energies achieved at RHIC and the LHC the produced 
QGP is a strong-coupled nuclear matter, which behaves like a relativistic fluid. 

2.1. Relativistic hydrodynamics and QGP viscosity 

In relativistic hydrodynamics (with no net conserved charge), one solves the con¬ 
servation equation for the energy-momentum tensor, 

= ( 1 ) 

where is given by 

T^''{x) = e[x)u^{x)vy{x) — P{x)A^'''{x) — Il{x)A^''{x) + 'k^''{x) . (2) 

Here e{x) and P{x) are the local energy density and pressure of the fluid. u{x) = 
7 [l,v(a;)] is the local four-velocity {u^u^ = 1), and A^''{x) = — u^{x)u'^{x). 

n(a:) and tt^'^{x) are the bulk pressure and the shear tensor. The first two terms 
are for ideal hydrodynamics, and the last two terms are viscous corrections. 
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Fig. 2. (Color online) Illustration a typical non-central heavy-ion collision in which the eccentric 
shape of initial collision geometry is translated into the elliptic flow anisotropy in the final state 
momentum space. 

One can see that there are 5 independent variables in ideal hydrodynamics (e, P 
and v), and 6 more unknowns in the viscous correction terms (11 and 5 independent 
components for The evolution equation for the energy-momentum tensor ren¬ 
ders 4 equations (for the energy density e and the fluid velocity v). The equation 
of state (EoS) will relate the pressure to the energy density, P = P{e), which is 
usually taken from lattice QCD calculation. Eor viscous hydrodynamics, additional 
evolution equations for the bulk pressure and the shear tensor are needed. They 
are usually obtained using the entropy principle, or derived from the kinetic the¬ 
ory via the moment method or the gradient expansion.l^SHlI] in the Israel-Stewart 
theoryPHHH we have: 




-h T:^''r]Td^ 

^ 'T— 


( 3 ) 


Here D = is the comoving derivative, -I- W'u^) — 

and V“ = A“^9;3 is the local spatial derivative. C, rj are bulk viscosity and shear 
viscosity, and rn, are the corresponding relaxation times. Currently most hydro¬ 
dynamics calculations ignore the bulk viscosity effects since shear viscosity provides 
the dominant contribution regarding the influence of the flow development in ultra- 
relativistic heavy-ion collisions.!^ 

In order to solve the relativistic hydrodynamics evolution equations, we need the 
spatial distribution of the energy momentum tensor at the starting time tq of the 
hydrodynamic evolution, which is usually obtained from various initial condition 
models. After solving hydrodynamic equations, one may obtain the spectra of the 
particles produced from the bulk matter according to the Coope-Erye description!^ 



( 4 ) 
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Fig. 3. (Color online) The p^ spectra for pions, kaons and protons in Pb-Pb collisions at 2.76 ATeV 
at the LHC. From top to bottom the curves are for 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 
50-60%, 60-70%, 70-80% collision centralities (for clearer presentation, the spectra are multiplied 
by a factor of 1000, 100, 10, 1, 0.1, 0.01, 0.001, lO”"^ and 10“^, respectively). Theoretical curves 
are from VISHNPE^and experimental data are from ALICE.I^The figures are taken from RefSU 


where is an infinitesimal element whose direction is normal to the hypersur¬ 

face E(a;), gi is the degeneracy factor, and fi{x,p) is the phase space distribution 
of particle species i. 

Fig. 1^ shows the pt spectra for identified hadrons (pions, kaons and protons) in 
Pb-Pb collisions at 2.76 ATeV at the LHC. Theoretical curves are from VISHNU 
calculation!^ and experimental data are from ALICE measurement s.^^ Here the 
QGP specific shear viscosity g/s in VISHNU calculations is taken to be a constant 
77 /s = 0.16 and the corresponding relaxation time is set as = 3g/{sT). The equa¬ 
tion of stateis taken from s95p-PCE parameterization which has been constructed 
by matching the lattice QCD calculation at higher temperature to a chemically 
frozen hadron resonance gas at lower temperature.^ The initial time for hydrody¬ 
namic evolution is set to be tq = 0.9 fm/c, and the initial entropy density profiles 
are generated from the MC-KLN model.^ One can see that except for the most 
peripheral collisions, where one does not expect hydrodynamics model to work well, 
the VISHNU model can provide a good description of the ALICE data for all three 
particle species over a variety of the collision centralities. 

Much of recent attention in the hydrodynamic studies of heavy-ion collisions has 
been paid to the quantitative extraction of transport coefficients such as the shear 
viscosity to entropy ratio g/s of the produced hot and dense QGP matter. Shear 
viscosity reflects the ability of the matter to transport momentum between different 
parts of the system. From kinetic theory, it is related to the mean free path A or 
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Fig. 4. (Color online) The integrated elliptic flow V 2 for all charged hadrons scaled by the eccen¬ 
tricity 62 in Au-Au collisions at 200 AGeV as a function of the charged hadron multiplicity per 
unit overlap transverse area (from RefPl). 


the transport cross section Ctr as: rj = |p(p)A = |(p)/crtr) where {p) is the average 
momentum carried by the matter constituents. Thus shear viscosity provides a 
direct measure of the interaction strength among the matter constituents. 

Fig.|^shows an effort to extract the specific shear viscosity rj/s from the charged 
hadron elliptic flow V 2 data in Au-Au collisions at RHIC.I^ The figure shows the 
Pt integrated elliptic flow V 2 scaled by the initial eccentricity €2 as a function of 
collision centrality [represented by the charged hadron multiplicity density per unit 
overlap area {1/S){dNch_/dy)]. The theoretical curves are from the VISHNU cal¬ 
culation^^ using two different initial condition models: MC-Glaubeil^ (left panel) 
and MC-KLINP (right panel). The experimental data show the elliptic flow mea¬ 
sured using two different methods: the event plane method (^ 2 ) and two-particle 
correlation methodu 2 { 2} l^lMl One can see that the centrality dependence of V 2 can 
be nicely reproduced by hydrodynamics calculation, and the presence of QGP shear 
viscosity suppresses the development of V 2 - In the figure, both theoretical curves 
and experimental data for V 2 are normalized by £2 to narrow down the sensitivity 
to the different experimental methods used for measuring V 2 - Since the magnitudes 
of £2 differ by about 20% in two initial condition models, the experimental data are 
shifted in two panels. Also due to the choice of different initial condition models, 
the extracted value of QGP rj/s has an uncertainty about a factor of 2-2.5. 

2.2. Initial state density and geometry fluctuations 

As has been shown, the knowledge of initial conditions, especially the geometry of 
the collision zone, is essential for the quantitative extraction of the shear viscosity 
of the QGP matter produced in heavy-ion collisions. For many years, most of the 
dynamical descriptions based on realistic hydrodynamics utilized averaged, smooth 
and symmetric initial density profiles (see Fig. [^. However, quantum fluctuations 
exist in the earliest stages of the collisions, such as the fluctuations of nucleon 
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Fig. 5. (Color online) The simulation of the fluctuating initial states in Au-Au collisions at 
200 AGeV at RHIC using a Glauber Monte-Carlo model.I^ These events are selected on pur¬ 
pose to obtain large values of eccentricities: £ 2 , £ 3 , £4 and £5 (from left to right). 


positions in nuclei, the fluctuation of color charges in nucleon, and so on. This leads 
to lumpy and asymmetric density profiles for the produced fireball which fluctuate 
one event to another, even for the collisions at a fixed impact parameter.l^^t^ 

To quantify the geometry (the geometric anisotropies) of the fluctuating initial 
states, for each event one may define the eccentricities, e„ = where e„ and 

are the magnitude and the orientation direction for the n-th order eccentricity 
vector. They may be calculated from the initial density prohles as follows: 




( 5 ) 


where p(r_L) is the (energy or entropy) density distribution in the transverse plane. 
Here the power index m is the radial weight, and often taken to be m = n for 
n > For n = 1 (dipole asymmetry), it is suggested to take m = 3.1^ One 

can see that for smooth initial conditions without fluctuations as shown in Fig. 
all odd Cn coefficients vanish. 

Fig. [5] shows a simulation of the fluctuating initial states produced in Au-Au 
collisions at 200 AGeV at RHIC using a Monte-Carlo Glauber model.l^ One can see 
that the density prohles of the initial state collision zones are lumpy and anisotropic 
in the transverse plane, and the initial conditions fluctuate violently from one event 
to another. The events shown in the hgure have been selected on purpose in or¬ 
der to obtain large values for £ 2 , £ 3 , £4 and £5 (from left to right), so the initial 
collision zones are mainly dominated by the elliptic, triangular, quadrangular, and 
pentagonal shapes, respectively. 

Fig. [ 6 ] shows the initial state £„ power spectrum for Pb-Pb collisions at 
2.76 ATeV at the LHC. Results are shown for four collision centralities, and three 
initial condition models (MC-Glauber,l^ MC-KL]>i^ and IP-GlasmeP^. One can 
see that in the ultra-central (0 — 0 . 2 %) collisions, all £„’s have roughly the same 
magnitudes since they all originate from the initial state huctuations. For less cen¬ 
tral collisions, since the collision zone has a pronounced elliptic shape, the values 
of £2 (and to some degree £ 4 , €q as well) are larger than the odd £„ coefficients. 
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Fig. 6. (Color online) The power spectrum of the initial state eccentricities e„ for Pb-Pb collisions 
at 2.76 ATeV. Results are shown for four collision centralities and three initial condition models 
(MC-Glauber, MC-KLN and IP-Glasma). The figure is taken from Ref.^ 


2.3. Hydrodynamic response and high-order anisotropic flow 

As a consequence of the interaction among the medium constituents, hydrodynamic 
evolution of the fireball will translate the initial geometric eccentricity into final 
state momentum anisotropy. To quantify such effect, one may perform the Fourier 
decomposition for the transverse momentum distribution of final state particles, 

— —oc 1 T 2 n{PT, y) COS - 4'n(PT, V))] , (6) 

where Vn and '^n are the magnitude and orientation angle (event plane) of the n-th 
order anisotropic flow vector v„ = Note that Vn and 'I'n are defined for 

a single collision, and can be pr (and y) dependent or integrated. The anisotropic 
flow may be obtained from the final state momentum distribution as follows: 

v„ = (7) 

Fig. [^show the scattered plots for {t 2 ,V 2 ) (left panel) and ( 53 ,^ 3 ) (right panel) 
from event-by-event hydrodynamics calculations for 20-30% central Au-Au colli¬ 
sions at 200 AGeV at RHIC.I^ One can see that the anisotropic flow coefficients V 2 
(U 3 ) show strong correlation to the initial state eccentricities €2 (ca)- This is true 
for their orientation angles as well: the final state event plane are strongly cor¬ 
related to the initial participant plane for n = 2 and n = 31^1^113 The linear 
response ^ 2/62 is stronger than v^/e^. For higher order harmonics, the linear corre¬ 
lation between and is spoiled by the non-linear contribution from lower-order 
harmonics . 1 ^ 
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Fig. 7. (Color online) Scattered plots for (€ 2 , 1 ^ 2 ) (left panel) and (e 3 ,i; 3 ) (right panel) from event- 
by-event hydrodynamics calculations for 20-30% central Au-Au collisions at 200 AGeV at RHIC 
(from Reffel). 
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Fig. 8. (Color online) The Vn power spectrum in ultra-central (0-0.1% centrality) Pb-Pb collisions 
at 2.76 AGeV at the LHcCH compared to viscous hydrodynamics calculations utilizing a few 
different initial condition models (from Ref. Ell). 
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Fig. § show the final state Vn power spectrum for ultra-central (0-0.1% central¬ 
ity) Pb-Pb collisions at 2.76 AGeV at the LHC compared to the viscous hydrody¬ 
namics calculations using a few different initial condition models One can see 
that in hydrodynamics calculations, higher order anisotropic flow harmonics are 
typically suppressed compared to lower order harmonics. This is true even in ideal 
hydrodynamics. We remind that all the initial eccentricities have roughly the same 
magnitudes in ultra-central collisions as shown in Fig. This is understandable 
since hydrodynamic evolution is more sensitive to larger scale structure (lower mo¬ 
mentum modes) of the initial density profiles, and less sensitive to smaller scale 
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Fig. 9. (Color online) The anisotropic flow coefficients Vn as a function of in Au-Au collisions 
at 200 AGeV at RHIC (left panel) and in Pb-Pb collisions at 2.76 AGeV at the LHC (right panel) 
(from RefPl), 


structure (thus higher order harmonics are suppressed). 

Fig. [9] shows the anisotropic flow coefficients as a function of px in Au-Au 
collisions at 200 AGeV at RHIC (left panel) and in Pb-Pb collisions at 2.76 AGeV at 
the LHC (right panel). The theoretical results are from (3-|-l)-dimensional viscous 
hydrodynamics calculations using the IP-Glasma-|-MUSIC model.l^ The data at 
RHIC are from PHENIX!^ and STAR.I^ The data points at the LHC are from 
ATLAS.I^The best descriptions to the experimental data give the average value of 
rj/s to be 0.12 at RHIC and 0.2 at the LHC. This means that on average the QGP 
medium produced at the LHC is less strongly coupled than that at RHIC. Since the 
temperature of the medium is higher at the LHC, this suggests that there is a strong 
temperature dependence for ry/s. The precise determination of the temperature- 
dependent specific shear viscosity rj/s(T) is one of the essential tasks in the current 
study of heavy-ion collisions. We note that with the use of the parameterization 
of temperature-dependent ri/s{T) from Ref.,l^a reasonable description of the flow 
data was also obtained from the IP-Glasma-|-MUSIC model.l^ 


2.4. Flow fluctuations and correlations 


Due to large fluctuations in the initial states, both e„ and Vn vary strongly from 
one event to another, even for a very narrow centrality bin. Since the averaged 
mainly reflect the hydrodynamic response to the averaged initial collision geome¬ 
try of the produced fireball, the measurements of event-by-vent distributions is 
essential in order to obtain direct insights into the fluctuations in the initial states. 
Both ATLAS and ALICE Collaborations have measured the event-by-event dis¬ 
tributions for Pb-Pb collisions at 2.76 ATeV at the LHCEHIMl Fig. 10 shows the 
results from ATLAS in several centrality bins. One feature of distributions is 
that they become broader from more central to more peripheral collisions due to 
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Fig. 10. (Color online) The event-by-event V 2 (left panel), D 3 (middle panel) and 1^4 (right panel) 
distributions for charged particles with px > 0.5 GeV in several centrality classes for Pb-Pb 
collisions at 2.76 TeV at the LHC, measured by ATLAs!^ 



Fig. 11. (Color online) Correlations between two different event planes (co^ 
j is an integer, and k is the least common multiple of m and n (from Ref.l^ 




"tn)]), where 


the increasing of the magnitudes. A strong centrality dependence is observed 
for the shape of V 2 distribution, while higher-order harmonics shows much smaller 
centrality dependence. The event-by-event distributions have also been studied 
in some details by a few groups utilizing relativistic hydrodynamics with various 
initial condition modelsE3[131Zl 

Various hydrodynamics calculations have shown that the elliptic flow V 2 and the 
triangular flow U 3 are mainly driven by initial state £2 and £ 3 . However, higher-order 
flow harmonics Vn {n > 3) arise from a combination of £„ and the non-linear mixing 
with lower-order harmonics. Such non-linear mode-mixing will manifest in the final 
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Fig. 12. (Color online) The vz{v 2 ) (top panels) and i’ 4 (i^ 2 ) (bottom panels) correlations measured 
in 0.5 < pt < 2 GeV in three centrality bins. In each panel, the correlation data are fit to functions 
that include both linear and non-linear contributions. The correlation data are also compared with 
re-scaled en(e 2 ) correlation from the MC Glauber and MC-KLN models in the same centrality 
bins. The figures are taken from Ref.1111 


state correlations between different order flow harmonics. Fig. [m shows the corre¬ 
lations between two different event-plane angles as a function of collision centrality. 
Experimental data are from Ref.^^ and theoretical curves are from Ref.^ with 
the use of two different initial condition models, MC-Glaubeil^ and MC-KLhP3 
(for later hydrodynamics evolution, ry/s is taken to be 0.08 and 0.2 for these two 
models, respectively). One can see that hydrodynamic calculation can qualitatively 
describe the observed centrality dependence of event plane correlations. We note 
that several final state event plane correlators exhibit very different centrality de¬ 
pendence from the initial state participant plane correlatorsThis is mainly 
due to the development of the mode-mixing between different flow harmonics as a 
result of non-linear hydrodynamics evolution of the fireball. 

Besides the correlations between different flow angles the magnitudes of 
different flow coefficients correlate to each other as well. Fig.[^shows the ATLAS 
measurement of the correlations between the magnitudes of different order flow 
harmonics {V3-V2 and V4-V2 correlations) using the event-shape selection method.!^ 
One can see that V 3 is anti-correlated to V 2 , which is very similar to the correlation 
between 63 and €2 in the initial states. This indicates V 2 -V 3 correlation reflect mostly 
initial geometry effects; this is expected since V 2 and V 3 are dominated by the 
linear hydrodynamics response. The V 2 -V 3 correlation data are fitted with a linear 
function V 3 = kv 2 + Ug, one can see that the magnitude of |A:|, which characterizes 
the strength of the anti-correlation, increases from central to peripheral collisions. 
This anti-correlation is understandable since when the overlap region becomes more 
elliptic, the fluctuation to large triangularity is constrained. For V 2 -V 4 correlation, 
both linear and non-linear collective dynamics contribute. The correlation data are 
well described by two-parameter fits with the form: v\ = + (ciu|)^, where the 
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Fig. 13. (Color online) The longitudinal correlation functions C2{^v) (1®^ panel) and CjjA r]) 
(right panel) in Pb-Pb collisions at 2.76 AGeV for three different centrality bins (from Ref.li^^l). 


linear term = cq is associated with £4 and the non-linear contribution = cit;| 
reflects the V 2 -V 4 mixing. The linear term depends weakly on centrality, whereas the 
non-linear term increases as the collisions become more peripheral^ (not shown). 


2.5. Longitudinal fluctuations 


Whe studying the anisotropic flow harmonics and final state correlations, one often 
utilizes a large pseudo-rapidity gap between correlated particles in order to minimize 
the contribution from non-flow effects such as resonance decays and jets. The use 
of large rapidity gap is reasonable when the initial density distribution and final 
state flow harmonics at different rapidities are prefectly correlated. However, the 
density profiles of initial states fluctuate not only in the transverse plane, but also 
in the longitudinal direction. Note that the fluctuations due to finite multiplicity 
for a given event are usually corrected using the sub-event methodEHIHIl Initial 
state longitudinal fluctuations may lead to the fluctuations and decorrelations of 
the final flow orientations at different pseudo-rapidities .IMtlooi One consequence is 
the reduction of the values of the final state anisotropic flows.^ 

To study the effect of longitudinal fluctuations, one may define the following 
correlation function between two rapidity bins A and HpSSl 


Cn{A,B) 


(Q„(A)Q*(B)) 

^(Q„(A)Q*JA))^(Q„(B)Q*JB)) ’ 


( 8 ) 


where Q„ = = (e™'^). In the continuum limit, Q is identical to flow vector 

v„. Fig. 13 shows the correlation functions C 2 (left panel) and C 3 (right panel) from 
relativistic hydrodynamics calculation for three different centrality bins for Pb-Pb 
collisions at the LHC.One can see that the anisotropic flows at different pseudo- 
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Fig. 14. (Color online) The Fourier coefficient V2{2} for hadrons, pions, kaons and protons as a 
function of py measured by ALfCE Collaboration using two-particle correlations in the 0 -20% 
multiplicity class after the subtraction of that from the 60-100% multiplicity class (from Ref. L^^ I). 


rapidities are not perfectly correlated, and the degree of such decorrelation becomes 
stronger with increasing pseudo-rapidity gap. The correlation C 2 of the second-order 
anisotropic flow shows a strong centrality dependence whereas the correlation C 3 of 
the third-order anisotropic flow is independent of centrality. This can be understood 
since U 3 arises purely from initial state fluctuations which is almost independent 
of collision geometry, while V 2 has a strong dependence on the collision centrality 
(the geometry of the collision zone). It was further shown that the longitudinal 
decorrelatons in the anisotropic flows are caused by the longitudinal fluctuations in 
the initial state density distributions with a twist structure as well as as additional 
random fluctuations on top of a twist.l^^ 

2.6. Anisotropic flow in small collision systems 

The relativistic proton-proton (p-p), proton-nucleus (p-A) and deuteron-nucleus 
(d-A) collisions are expected to provide the baseline for studying the transport 
properties of the hot and dense QCD matter produced in A-A collisions. The col¬ 
lective behavior and anisotropic flow are expected to be much weaker in these 
colliding systems since the size of the produced matter is much smaller (even 
though a mini quark-gluon plasma is produced). However, ALICE, ATLAS and 
CMS Collaborations have observed a clear collective behavior in p-Pb collisions at 
5.02 ATeV,E2iHin3| very similar to peripheral Pb-Pb collisions at 2.76 ATeV at the 
LHC. Similar correlation results have been obtained for d-Au collision at RHICPSl 
and high multiplicity p-p events at the LHC energies 

The origins of the observed correlations and collective behavior in these smaller 
colliding systems are still in debate. One natural explanation for the observed col- 
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lectivity and correlations is the hydrodynamic expansion of the fireball with fluctu- 
ating initial conditions jlflZlIIini which suggests that a mini-QGP may be produced 
in these small colliding systems. A strong support for hydrodynamics explanation 
is the mass ordering of the elliptic flow V 2 observed by the ALICEJi^ (see Fig. 
[l^ , which was reproduced later on by hydrodynamic calculationsEnilllll Recently 
Refllisi shown that the long-range azimuthal correlations in both p-p and p- 
Pb collisions can be explained by the incoherent scattering of partons using the 
AMPT model with string melting.^^Sl observed two-particle azimuthal corre¬ 
lations was also fitted by the two-gluon emission mechanism in the initial states 
within the framework of Color Glass Condensate.^^^^ Recently relativistic ^He-Au 
collisions were proposed in order to exploit the intrinsic triangular geometry in the 
initial satesEi^lii^ By detailed and systematic comparisons of relativistic p-A, d- 
A, ^He-A and A-A collisions, it is promising that one can disentangle the initial 
geometry contribution from the final state viscosity effect. 


3. Color Opacity of QGP 


Hard partonic jets that are produced from early stage scatterings provide important 
probes to study the properties of the QGP created in high energy nuclear collisions. 
During their propagation through the hot and dense nuclear matter, they interact 
with the medium constituents via elastic and inelastic collisions, and usually lose 
energy in the processThis is often referred to as jet quenching. The main 
purpose of jet quenching studies in heavy-ion collisions is to understand the detailed 
mechanisms of jet-medium interaction from which one may infer useful information 
about the produced hot and dense QGP. 

One of the important jet quenching observables is the suppression of the single 
inclusive hadron yield at high transverse {pt) in A-A collisions as compared to that 
from elementary p-p collisions ]122M221 Xo quantify such effect, one may define the 
nuclear modification factor Raa as follows: 


RAA{PT,y,'<l^) 


1 dNAA/dpTdydip 
Ncoii dNAA/dprdyd-ijj ’ 


( 9 ) 


where iVcoii is the average number of binary collisions for a given collision centrality 
bin. Fig. [l5| shows the nuclear modification factors Raa for charged hadrons, direct 
photons and other particles; the left panel shows PHENIX measurements for most 
central Au-Au collisions at 200 AGeV at RHIC, and the right shows CMS measure¬ 
ment for most central Pb-Pb collisions at 2.76 ATeV at the LHC. One can see that 
high Pt hadron yields are strongly suppressed in A-A collisions compared to binary 
collision scaled p-p collisions. In contrast, Raa for high pt photons is consistent 
with unity. This means that the observed strong nuclear modification for high pt 
hadron production is a final state effect from the interaction of hard partonic jets 
with the hot and dense medium. 

There are many other experimental signatures (observables) for jet quenching 
at RHIC and the LHC, such as the modification of the correlated back-to-back 
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Fig. 15. (Color online) Left: Raa for hadrons, direct photons and no n-ph otonic electrons in most 
central Au-Au collisions at 200 AGeV at RHIC, measured by PHENIX.li^ Right: Raa for hadrons, 
direct photons, a nd Z /W bosons in most central Pb-Pb collisions at 2.76 ATeV at the LHC, 
measured by CMS.l^^ 


dihadrons and high px photon-hadron pairs,the suppression of reconstructed 
jet and dijet production,E2nHIl3l the nuclear effect on jet fragmentation and jet 
shape functions liMtlMI Similar to the single inclusive hadron Raa, one may define 
the nuclear modification factors for these jet quenching obervables. 

3.1. General framework for jet quenching study 

In perturbative QCD, processes that involve large momentum transfer can be de¬ 
scribed as the convolution of the parton distribution functions (PDFs), hard par- 
tonic scattering process, and the final state fragmentation function (FFs). For ex¬ 
ample, the cross section for the single inclusive hadrons at high px in p-p collisions 
may be calculated as followsfi^ 

dapp^hX ~ ^ ^ dXa dXf) j dzj fa/pi^Xa, gf') 0 fb/p{^b, 

abj 

(g) doab^jxigf, gp, gn) ^ gp), (10) 

where fa/pi^a, gf) and fb/pixb, gf) are two PDFs with Xa and Xb the momentum 
fractions of the incoming partons, doab^jx is the parton scattering cross section, 
and Dj^h{zj, gp) is the FF for the parton j to the hadron /i, with Zj is the momen¬ 
tum fraction of the outgoing hadron. There are three momentum scales involved 
here: the factorization scales gf and gp and the renormalization scale gp, they are 
usually taken to be the same as a typical hard scale (Q) involved in the process, 
such as the hadron pp. The PDFs and FFs are non-perturbative and universal func¬ 
tions, and obey the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations 
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(GeV/c) 


Fig. 16. (Color online) The cross section for inclusive hadrons at high p'r in p-p collis ions at RHIC 
measured by PHENIX, compared to NLO perturbative QCD calculations (from Ref.^^. 


for their scale evolutionsThey are usually determined by the global fit to 
e+e” experiments, deep inelastic scatterings (DIS) and p-p collisions, etc. 

Fig. shows the production cross section for single inclusive hadrons at high 
Pt in elementary p-p collisions at RHIC energies. One can see that the high pT 
hadron production can be described quite well by the next-to-leading order (NLO) 
pertubative QCD calculations.*!^ This indicates the properties of hard jets in vac¬ 
uum are well understood from both experimental and theoretical sides, which serves 
as the baseline for studying jet-medium interaction and the nuclear modification of 
hard jets in relativistic heavy-ion collisions. 

When studying jet quenching in relativistic heavy-ion collisions, two nuclear 
effects need to be taken into account. One is the nuclear modification of PDFs, 
i.e., the PDF in nucleus fa/A is different from the free proton PDF fa/p', this is the 
initial state nuclear effect, and often called cold nuclear matter (CNM) effect. CNM 
effect can be included by defining the nuclear modification factor Rf{x, Q^) for the 
PDF: R^(x,Q^) = fa/Aix,Q^)lfalpix,Q^), which is often obtained by the global 
fit to DIS, p-A and d-A collisions. Currently several parameterizations of PDF 
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Fig. 17. (Color online) Typical diagrams for collisional (left panel) and radiative (right panel) 
energy losses of a hard parton propagating through a nuclear medium. 


nuclear modification factors Rf{x,Q'^) are available (e.g, EPSOOj^^ HKNO?,*!^ 
and 

The second effect is due to the production of the hot and dense QGP medium in 
relativistic heavy-ion collisions, which we may call hot nuclear matter (HNM) effect. 
Partonic jets produced from the initial hard scatterings have to travel through and 
interact with the produced QGP before fragmenting into final observed hadrons. 
With the inclusion of both cold and hot nuclear effects, the single inclusive hadron 
production in heavy-ion collisions may be calculated as follows: 

dUAB^hX « X! fa/A(.Xa) G fb/Bi^b) G dUab^jX ® ® ,(11) 

abjj' 

Here {pj' \pj) describes the hot nuclear matter effect, i.e., the interaction of the 
hard partons j with the colored medium. It is noted that although the above fac¬ 
torized formula has been widely used in phenomenological studies of jet quenching 
in heavy-ion collisions, there is currently no formal proof of such factorization yet. 
The assumption of factorization is consistent with various jet quenching studies. 

3.2. Radiative and collisional jet energy loss 

Hard partons may lose energy in hot and dense nuclear medium via a combination of 
elastic collisions with the medium constituents and inelastic radiative process. Fig. 
|17| shows two typical diagrams for calculating parton energy loss originating from 
elastic (left panel) and radiative (right panel) processes. The energy loss experienced 
by the primary parton in 2 —> 2 elastic collisions with the medium constituents is 
usually called collisional or elastic energy loss. The additional in-medium radiation 
induced by multiple scatterings may take away a fraction of energy from the primary 
parton, which is usually called radiative energy loss. 

Radiative energy loss process has been regarded as the most important compo¬ 
nent in studying jet quenching in relativistic heavy-ion collisions. Even in vacuum, 
hard partons produced from the early stage scatterings will undergo splitting pro¬ 
cesses. In a hot and dense nuclear medium, the parton splitting processes will be 
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modified due to the rescatterings of the propagating parton with the medium con- 
stituents. A number of parton energy loss approaches have been developed to study 
medium-induced radiative process, namely Baier-Dokshitzer-Mueller-Peigne-Schiff- 
Zakharov (BDMPS-Z) lliMUZl Gyulassy, Leva! and Vitev (GLV),li^2Hl50l Amesto- 
Salgado-Wiedemann (ASW) higher twist (HTjPSSHHH Arnold-Moore- 

Yaffe (AMY^PSfiHUH formalisms. One may refer to Ref.l^^for a detailed comparison 
of different parton energy loss formalisms. 

In the study of radiative parton energy loss, one key quantity is the single emis¬ 
sion kernel (e.g., the differential gluon radiation spectrum dNg/duidk\dt). Various 
approaches make different assumptions about the traversed medium and utilize dif¬ 
ferent methods for treating multiple scatterings. One important effect in medium- 
induced emission or splitting processes is the Landau-Pomeranchuk-Migidal (LPM) 
effect For collinear or small angle radiation, a finite period of time is re¬ 

quired to complete the radiation process; this time is called the formation time 
Tf ^ 2a;/fc^, with u) and fc_L the energy and transverse momentum of the radiation. 
If the formation time is larger than the mean free path A of the propagating parton, 
the multiple scatterings on the propagating parton can no longer be treated as in¬ 
dependent. Such quantum interference between successive scatterings is called the 
LPM effect which will lead to the suppression of the radiation spectrum compared 
to the Bethe-Heitler incoherent multiple scattering limit. For the QCD case, since 
the radiated gluons carry color charge, the medium modification of the radiation 
spectrum is more dominated by the rescatterings on the emitted gluons. 

To obtain the multiple gluon emission, one common practice is the repeated 
application of single gluon emission kernel. Such recipe neglects the interference 
between different emissions which should be included in a full calculation of multiple 
parton final state. One popular method is the Poisson ansatz which has been widely 
used in GLV and ASW calculations.li^^llllll xhe key quantity is the probability 
distribution P{AE) of parton energy loss, which may be obtained as follows: 


P{AE) = £ 


.-{Ng) 


n —0 


n / 


djjj 


5 ( AE — iOi 


( 12 ) 


where dN/duj is the single gluon emission spectrum and (Ng) = J dujdNg/du! is the 
mean number of radiated gluons. In the AMY formalism, the following rate equa¬ 
tions are solved for the parton momentum distributions f{p) = dN (p) / 

dt J 


f{p + k,t) 


dT{p -h A:, A:, t) 

dkdt 




dT{p, A:, t) 
dkdt 


(13) 


where dT(p, k, t) j dkdt is the transition rate for a parton with momentum p to lose 
momentum k. In the HT formalism, one solves the following DGLAP-like evolution 
equations for the medium-modified fragmentation function D{z, Q^)!iSSltiIU 


dD{z,Q‘^,q ) as f dy 


9 In <5^ 


27T 


P{y) / dC K{C ,Q\q ,y)D{z/y,Q\q y), (14) 
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Fig. 18. (Color online) Right: Average energy loss of charm (left panel) and bottom (right panel) 
quarks passing through the nucl ear m edium created in most 0-7.5% central Pb-Pb collisions at 
2.76 ATeV at the LHC (from Ref.lHl), 


where q is the jet light-cone energy, y the energy fraction of the radiation, and 


is the jet location. P{y) is the vacuum splitting function and K{(~,y) 
is the parton-medium interaction kernel. 

Collisional energy loss of hard partons propagating through the nuclear matter 
was hrst studied by Bjorken.^^ Compared to medium-induced radiation, collisional 
energy loss is usually considered to be small for light flavor (leading) partons, es¬ 
pecially at high energies imiElSEZll However, collisional energy loss may give siz¬ 
able contribution to nuclear modification factor Raa at RHIC and the LHC en- 
ergiesEUllIlZI In contrast, elastic collisions are usually considered as the dominant 
mechanism for heavy quark energy loss, especially at low and intermediate energy 
regimes liZilliZSl This is because the phase space for collinear medium-induced radi¬ 
ation is reduced by the presence of the finite masses of heavy quarks (usually called 
the dead-cone effect?^. This mass effect will diminish when going to high energy 
regimes where heavy quarks become ultra-relativistic as well and behave more like 
light flavor partonsThis can be seen in Fig. 18 which compares the radiative 
and collisional contributions to the energy loss of charm (left panel) and bottom 
(right panel) quarks. We note that collisional energy loss may play an essential role 
in the study of the nuclear modiflcation of full jetJ^^ and the medium response to 
the propagation of hard jetsjUUEll these will be discussed in later subsections. 


3.3. Jet quenching phenomenology at RHIC and the LHC and 
quantitative extraction of q 

In recent years, various phenomenological studies on jet quenching have been per¬ 
formed for a wealth of experimental observables. Much effort has been made to 
the quantitative extraction of various jet transport parameters,^^ such as the 
transverse momentum diffusion rate q = d{Ap'^_^)/ dtpJ^ the elastic energy loss rate 
e = dE/dt « dp\\/dt, etcESSESH At leading order (LO), these transport coefficients 
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Fig. 19. (Color online) The values of scaled jet transport parameter q/T^ extracted by JET 
Collaboration by using single inclusive hadron suppression factor Raa a-f both RHIC and LHC. 
The values are for a 10 GeV quark jet at the cen ter o f most central collisions at an initial proper 
time TQ = 0.6 fm/c. The figure is taken from Ref.USU 


quantify the transverse and longitudinal momentum transfers experienced by the 
propagating partons via 2 —> 2 elastic collisions with the medium constituents. 

Recently a significant collaborative effort was carried out within the framework 
of JET Collaboration; systematic phenomenological studies were performed on the 
experimental data on the nuclear modification of single inclusive hadrons at large 
Pt in heavy-ion collisions at the RHIC and the LHC.li^ Five different existing ap¬ 
proaches to parton propagation and medium-induced energy loss in dense medium 
were used for this systematic survey: McGill-AMY jUIElll Martini-AMYHT- 
Mplo] ht-bw,e^ DGLV-CUJET.E^ The space-time profile of the QGP medium 
was simulated by the (2-|-l)-dimensional or (3-|-l)-dimensional hydrodynamic mod¬ 
els. The goal of this systematic study is to use the constraint from the experimental 
data and to quantitatively extract the jet quenching parameter q and its systematic 
uncertainties (including model dependence). 

By comparing each model calculation with the experimental data and hxing the 
model parameters, the effective jet transport coefficient q may be obtained. Fig. |19| 
shows the extracted values of q scaled by for a quark jet with energy 10 GeV at 
the highest temperatures reached in the most central Au-Au collisions 200 AGeV 
at RHIC and Pb-Pb collisions 2.76 ATeV at the LHC. Using the measured single 
hadron Raa, the range of the values for q/T^ is obtained as: 

q _ r3.5±0.9, Tw 370 MeV (at RHIC), 

CsT^ “ \ 2.8 ± 1.1, T « 470 MeV (at the LHC). 


(15) 
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This gives q « 1.2 GeV^/fm at RHIC and q « 1.9 GeV^/ fm at the LHC. 

One can clearly see the strong temperature dependence for the scaled parameter 
q/T^, which may be explored by extending the current study to future higher energy 
Pb-Pb collisions at the LHC and lower energy collisions at RHIC. The expected 
values of q/T^ at 0.063 ATeV, 0.130 ATeV and 5.5 ATeV are shown in the figure. 
Also shown is the value of (Jn/T^s nuclei, extracted from jet quenching 

studies in DIS.l^^ One can see that the values of jet parameter q in hot QGP 
medium are much higher than those in cold nuclei. The values of q/T^ from a NLO 
AdS/CFT calculation is also shown for comparison and one can see that they 
are within the range of q values from JET Collaboration. Note that the SYM values 
quoted here have been obtained for asvM = 0.22-0.31 and included the effect due 
to different numbers of degrees of freedom in SYM and QCD. 

In the future, one may include more theoretical model calculations and utilize 
more jet quenching observables in such systematic phenomenological studies.li^ 
Besides q, other jet transport coefficients (such as e) may also play significant roles 
in jet-medium interaction and medium-induced parton energy lossE^lilSiH^lli^ 
To estimate the systematic uncertainties in the LO calculations, it is essential to 
develop a fully NLO framework for studying jet propagation and modification in 
dense nuclear matter. Recently the NLO radiative correction to transverse momen¬ 
tum broadening and the renormalization of jet quenching parameter q have been 
investigated in ftef.llS^lOnil ^11 these important ingredients should be included in 
the future study of jet quenching in heavy-ion collisions. 

3.4. Jet quenching from Lattice QCD 

While lattice QCD cannot provide a full description of the dynamical evolution of 
heavy-ion collisions, it can provide important guidance for phenomenological jet 
quenching studies, e.g., many transport coefficients (such as q) may in principle be 
computed using lattice QCD. Ref .*2221 carried out the first lattice QCD calculation of 
jet quenching parameter q using the operator product expansion method; the phys¬ 
ical q is related to an infinite series local operators in an unphyical regime of jet mo¬ 
menta via dispersion relations. The calculation was performed in quenched SU (2) 
and the extension to 517(3) with 2 flavors of quarks produced g=I.3-3.3 GeV^/fm 
for a gluon jet at a temperature T = 400 MeV, which is similar to the values from 
JET Collaboration. Recently within a dimensionally reduced effective theory (elec¬ 
trostatic QCD), Ref.lSMl obtained a value of (j=6 GeV^/fm at RHIC energies, which 
is about twice the value compared to NLO perturbative QCD calculation.l^SSI 

We note that the jet quenching parameter q is the second moment of the prob¬ 
ability distribution P{kj_,L) for the transverse momentum transfer, 



(16) 


More detailed information about jet-medium interaction are contained in the full 
distribution of the probability distribution P{k±) or the transverse collision kernel 
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Fig. 20. (Color online) The collision kernel C{kj_) (left panel) and the probabilit y di stribution 
P(fcx, L) (right panel) compared with the LO and NLO calculation results from ~ g^T 

is the effective coupling of the dimensionally reduced effective theory and 0q is a dime nsion less 
number related to the lattice spacing: = 2Nc/{g‘^Ta). The figures are taken from RefP^ 


C{k±). In Ref. ,12221 a first-principle calculation of the collision kernel C{k±) was 
carried out (see Fig. [20| ). At small transverse momentum kj_, the shape of C{kj_) 
was found to be consistent with Gaussian distribution. Since the calculation only 
includes the soft modes, it breaks down when k± becomes large. 

The Gaussian form for P{k±) and C{k±) at small k± is expected from both 
weakly-coupled and strongly-coupled calculations. However the large k± behaviors 
are very different in two scenarios: for weakly coupled QGP, the power tail for 
large k± is proportional to 1/fc^ , whereas in Af = 4 SYM theory (in the large- 
Nc and strong coupling limits), there is no power tail at all .12211 One interesting 
task is whether the full distribution of the collision kernel can be obtained from 
phenomenological jet quenching studies (as well as from lattice QCD calculation). 

3.5. Full jet evolution and energy loss 

The basic idea of full jet reconstruction is to recombine final state jet fragments and 
obtain the information about the original hard partons and investigate the medium 
modification effects. Full jets are expected to provide more differential information 
about jet-medium interaction than leading hadron observables since they include 
both leading and sub-leading jet fragments. Early full jet studies in heavy-ion col¬ 
lisions were performed in Au-Au and Cu-Cu collisions at RHIG by both STAR 
and PHENIX CollaborationsIn spite of large experimental uncertainties 
and strong dependence on jet reconstruction algorithms, substantial nuclear mod- 
ihcation of full jets in heavy-ion collisions have been observed as compared to p-p 
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Fig. 21. (Color online) Left: The energy asymmetry Aj = (Et,i ~ ^T, 2 )/(®T.l + ®T, 2 ) distribu¬ 
tion for dijets in p-p and Pb-Pb collisions at 2760 AGeV at the LHC (from Ref.liSU), Right: The 
energy asymmetry zj-^ = Ej’ jjEx^-y distribution for corr elate d isolated photons and jets in p-p 
and Pb-Pb collisions at 2760 AGeV at the LHG (from Ref.l^iiJ) . 


collisions. The strong dependence on the jet resolution parameter R for the nu¬ 
clear modification of jet production indicated the broadening of full jets due to the 
interaction between jets and the hot and dense nuclear medium .12121 

The launch of the LHC has increased the center of mass energy by more than 
a factor of 10 compared to the top collision energies at RHIC, which enables us 
to investigate the medium modification on the propagation of jets with transverse 
energies over a hundred GeV. The first measurements of the nuclear modification of 
full jets at the LHC is the correlated back-to-back jet pairs. We observed a strong 
centrality dependence for the modification on dijet transverse energy imbalance 
Aj = {Ets — Et, 2 )/{Eta + Et, 2 ) distribution in Pb-Pb collisions as compared to 
p-p collisions, while the distribution of their relative azimuthal angles is largely 
unmodified. Similar results have been obtained for the full jets correlated with high 
Pt direct photons in Pb-Pb collisions at the LHC.*^^ These results indicate that 
the subleading jets may experience a significant amount of energy loss after passing 
through the produced hot and dense QCD matter. Various jet energy loss calcu¬ 
lations have been performed to explain the nuclear modification of the transverse 
energy imbalance between correlated dijets and photon-jet pair^HEUHlIZI (e.g. see 
Fig.p). 

Fig. shows a schematic illustration of the evolution of a full jet in a QGP, 
where the thick solid arrowed line through the center represents the leading parton, 
and other lines represent the accompanying radiated gluons. Compared to leading 
harden observables, we need to consider a few additional ingredients when studying 
full jet evolution and energy loss. Besides the primary parton, the radiated gluons 
may interact with the medium and lose energy in the process. The radiated gluons 
may get deflected by the medium constituents as well; some of the radiations may 
be scattered out of the jet cone. Therefore, the total energy loss of the full jet is 
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Fig. 22. (Color online) Illustration of the evolution of a full jet in a quark-gluon plasma. Various 
medium-induced processes that contribute to full jet energy loss and medium modification are 
shown. 


the sum of elastic energy loss experienced by the leading parton and the radiated 
gluons, together with the gluons that are kicked out of the jet cone. 

In the spirit of the above picture, Ref.*^^ performed the first quantitative anal¬ 
ysis of physical processes that are responsible for full jet energy loss. It solves the 
following transport equation for the radiated gluons: 


dfg{uj,k±,t) 

dt 





dujdk^dt ’ 


(17) 


where fg{uj,k±,t) is the three-dimensional momentum distribution of the accom¬ 
panying gluons. The first and second terms in the above equation describe the 
evolution of the radiative gluons which may transfer energy into the medium via 
elastic collisions and accumulate transverse momentum in the process. The last 
term is a source term which represents the contribution from the medium-induced 
gluon radiation. After they are produced, the medium-induced radiative gluons 
will interact with the medium, lose energy and accumulate transverse momentum 
during their propagation. 

By solving the above equation, one may obtain the information about the full 
jets after passing through the medium and calculate the energy loss from the jet 
cone. One may decompose the energy of the original full jet into a few parts: 


Ajet — E/in “t” E/iost — E/in,rad T E/out,rad T E/out,brd 4” E/t}i coll ■ (1^) 

It is interesting that the most significant contribution is found to be the collisional 
energy loss experienced by the radiated gluons Elth.coiif^^ The transverse momen¬ 
tum broadening of the radiative gluons Elout.brd also gives a sizable contribution 
to full jet energy loss. Since the radiation is mainly dominated by small angle ra¬ 
diation, the contribution from the radiation directly outside the jet cone Elout.rad 
is small. This finding is not surprising since one expects the soft components of 
the jet or the accompanying gluons at large angles to experience stronger medium 
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modification than the inner hard core of the jet. Such effect is also referred to as 

jet collimation.l21SI 

Similar picture for full jet energy loss has been obtained in RefElHIElS It is ar¬ 
gued that if the medium color field A varies over the jet transverse size the 
color coherence of the shower partons may be destroyed by the interaction with 
the medium. Such color decoherence effect may greatly increase the phase space for 
soft and large angle radiation, as compared to the traditional BDMPS-Z radiative 
energy loss formalism. Treating multiple gluon emissions as a probabilistic branch¬ 
ing process, one may solve the following rate equation for the gluon momentum 
distribution D{x,t) = xdNg/dx, 


dD{x, t) 

57 




r) 


z 



D{x,t) 


(19) 


Two terms in the above equation are gain and loss terms. Focusing on the small-a; 
behavior, one may take k{z) = 1. Given the initial condition D(x, r = 0) = i5(a: — 1), 
one may obtain the solution for the gluon distribution: D{x,t) « 

total energy contained in the spectrum may be calculated as, e(r) = dxD{x, t) « 
, and one can see that it is not conserved. Therefore with increasing time, a 
substantial fraction of jet energy could be lost from the gluon spectrum. 

To include the contribution from radiation at large angles, one may introduce a 
scale Xq, below which the radiation is directly outside jet cone. One may further in¬ 
troduce a scale Xth = F !E to include the contribution that when the gluon energies 
are as low as as the medium energy scale, the radiated gluons may “thermalize” 
and disappear from the jet. So the total phase space for the branching gluons (and 
the total energy of the original jet) may be divided into three partsP^ 


Ejet — Ejn “t“ Elost — E{n^{x ^ :ro) -t- Eout(3:'th ^ x ^ :ro) -1- Efiow(^ ^ ^th)- (^0) 


The branching gluons with x > xq are inside the jet cone. For Xth < x < xq, the 
radiations are directly outside the jet cone. For x < a;th, the gluons will thermalize 
and flow into the medium. The largest contribution to full jet energy loss is found 
to be the energy flowing out of the gluon spectrum. The radiation directly outside 
jet cone gives small contribution. 


3 . 6 . Full jet substructure 


To fully characterize the jet quenching effects, we need to study not only the total 
energy of the full jets, but also the modification of the jet substructures. Jet shape 
function p(r) describes the radial distribution of the momentum carried by the jet 
fragments. It is defined as follows: 

1 1 
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( 21 ) 


where = ^(??/t — ?7jet)^ + i't’h — </'jet)^j and 6r is the bin size. The differential 
jet shape function is normalized to unity, J p{r)dr = 1. Jet fragmentation function 
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CMS, = 2.76 TeV pp, J L dt = 5.3 pb*' PbPb, J L dt = 150 pb 
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Fig. 23. (Color online) Top: The jet shape functions p(r) in Pb-Pb collisions at 2.76 ATeV for jets 
with > 100 GeV/c in five centrality bins. Bottom: The nu clear modification factors for jet 
shape functions. The figures are taken from CMS measurements.^^ 


D{z) provides the information about the momentum spectrum of the jet fragments. 
It is defined as follows: 


D{z) 


1 dNh 
iVjet dz 


( 22 ) 


where z = • Pt*/IPt*P momentum fraction of the jet fragments. 

Fig. 23 shows the nuclear modification of the differential jet shape function p{r) 
in Pb-Pb collisions at the LHC measured by CMS Collaborations.!^^ The jet shape 
functions in peripheral Pb-Pb collisions are similar to those in p-p collisions while 
in more central Pb-Pb collisions we observe an excess at large radius r > 0.2 and 
a depletion at intermediate radii (0.1 < r < 0.2). This indicates the broadening of 
the full jets after they pass through the medium. At very small radii, little change is 
observed for jet shape function, suggesting that the energy distribution in the inner 
hard core of the jet is not affected by the jet-medium interaction. These observations 
are consistent with the previous CMS finding that the lost energy from the jets is 
found at large distances from the jet axis.*^^ 

Fig. shows the jet fragmentation functions p-p and Pb-Pb collisions at the 
LHC measured by CMS Collaboration.!^^ Compared to the earlier CMS measure¬ 
ment, !i^the charged particles with lower values of pr have now been included in the 
full jet reconstruction. We observe a clear nuclear modification of the jet fragmen¬ 
tation function in Pb-Pb collisions, which grows with increasing collision centrality. 
In peripheral collisions (50-100% centrality bin), the ratio of PbPb/pp is almost 
flat at unity. For the most central collisions, we observe a significant excess at high 
^ = ln(l/z) (low z) and a depletion at intermediate These results indicate that 
the spectrum of the particles inside the full jets receives an enhancement in the soft 
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Fig. 24. (Color online) Top: The jet fragmentation functions D{z) in Pb-Pb collisions at 2.76 ATeV 
for jets with 100 < Pt < 120 GeV/c in four Pb-Pb centrality bins. Bottom: The nuclear modifi¬ 
cation factors for jet fragmentation functions. The figures are taken from CMS measurements!!^ 


regime, compared to p-p collisions. ATLAS have also measured the nuclear mod¬ 
ification of jet fragmentation function and obtained similar results: an enhanced 
yield of low and large z fragments together with a suppressed yield of intermediate 
z fragments. 

The above observations are qualitatively consistent with jet energy loss calcu- 
lationsEmHUHllll There exist several possible sources that could affect the nuclear 
modification of jet substructure and contribute to the redistribution of jet energy 
inside the cone. The broadening of the parton shower may be contributed from 
both elastic collisions with medium constituents and the medium-induced radia¬ 
tions. The induced radiation may lead to the softening of the momentum spectrum 
of the jet fragments as well. Particles produced from the medium response to jet 
transport may contribute to the redistribution of the energy inside the jets.l^LO Djf. 
ferent hadronization mechanisms may lead to different jet fragmentation profiles.!2^ 
Further detailed studies on the nuclear modifications in different r of jet shape dis¬ 
tribution and different z of jet fragmentation function should be able to put tighter 
constraint on the modeling of jet-medium interaction. 

3 . 7 . Medium response to jet transport 

Jets may lose energy via a combination of elastic collisions and inelastic radiative 
processes when propagating through a QGP medium. Some of the lost energy from 
the hard jets is deposited into the medium, and may contribute to many jet-related 
observables. Various studies have shown that the propagating hard jets may excite 
the medium and produce Mach cone structureIf such structure is observed 
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Fig. 25. (Color online) The medium response to the energy deposited by a p rimary hard quark 
(left panel) or by a quark-initiated parton shower (right panel) (from RefPSll). 


in relativistic nuclear experiments, it will provide a direct probe to the speed of 
sound of the produced nuclear matter. However in heavy-ion collisions, the Mach 
cone pattern may be distorted by the large collective flow developed during the 
hydrodynamic expansion of the bulk matter The cone structure is also 
sensitive to various transport properties such as the shear viscosity to entropy 
density ratio rj/s of the QGP mediumlUHElI] 

One approach for studying the medium response is the use of full Boltzmann 
transport codes in which the jet propagation and the medium response are simu¬ 
lated at the same time within a single Monte-Carlo packageIn this method, 
the bulk matter is modeled by a collection of quasi-classical partons, and the jet- 
medium interaction is treated the same way as the interaction among the medium 
constituents. Another approach is to solve the following hydrodynamic equations. 


= r{x) 


dE dp± dpi\ \ 
dt ^ dt ’ dt ) 


(23) 


Here represents the rate of the energy and momentum deposited by the 

propagating hard jets. The source term J^(x) at a given space-time location may 
be obtained from jet energy loss calculations, e.g., the energy deposited into the 
medium is equal to the collisional energy loss from the propagating jet. For a single 
parton in a high temperature QGP medium, the collisional energy loss (deposition) 
rate in the leading logarithmic approximation may be obtained as follows: 


dE 

dt 


dep 


dE 

dt 


coll 



iET\ 
ml ) 


(24) 


where ml = 47ras(I -I- Nf/6)T‘^ is the Debye screening mass squared. 

It should be noted that jets are composed of collimated showers of partons. 
Therefore, not only the leading parton, but also the radiative emissions can de- 
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Fig. 26. (Color online) The medium response generated by a parton shower originating from a 
primary hard quar k wit h the average gluon emission angle being 6 = 0.1 (left panel) and 9 = 0.7 
(right panel), (frorrPUl). 


posit energy and momentum into the medium via scattering with the medium 
constituents. In Ref. plHlElll it is found that the length dependence of energy depo¬ 
sition rate by a jet shower is significantly enhanced compared to that by a single 
parton. The medium responses to the energy deposited by a single parton (left 
panel) and by a parton shower (right panel) are compared in Fig. 25 The Mach 
cone structure is seen for both energy deposition cases, but the conical pattern is 
strongly enhanced for a parton shower. Another important factor is the spatial dis¬ 
tributions of the parton shower and the energy and momentum deposition profiles. 
Fig. compares the medium responses to a parton shower with small emission 
angle 0 = O.I (left panel) and with large emission angle 0 = 0.7 (right panel).!^^ 
One can see that for nearly collinear emissions, a nice Mach cone structure can be 
produced. In contrast, a well-defined Mach cone is hardly seen for the case of larger 
angle emissions; the medium response is now more like a superposition of several 
energy density perturbations. 

In the above two studies (Fig. 25 261, the medium is taken to be static, with 
a constant temperature (density). In realistic event-by-event simulations which in¬ 
clude both initial state and parton energy loss fluctuations, the energy and momen¬ 
tum deposition profiles may vary strongly from one event to another.^^MI Since the 
realistic medium probed by the propagating jet is dynamically evolving, expanding 
and cooling, the energy deposition rate will hrst increase and then decrease as a 
function of evolution time. We also note that in many model calculations, a cutoff 
energy is often used to determine which part of radiation phase space “thermalize” 
into the medium. The change of such cutoff (separation) scale may lead to quite 
different energy loss and deposition profiles for the propagating jetlHolEll 
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4. Summary 

One of the main goals of high energy nuclear collisions is to create the hot and 
dense nuclear matter with (energy) densities well above the normal nuclear mat¬ 
ter and study its various novel properties. In this report, we have provided some 
basic information and recent progresses on the study of quark-gluon plasma using 
relativistic heavy-ion collisions at RHIC and the LHC. Our focus was given to the 
anisotropic collective flow and jet quenching phenomena, which are two most im¬ 
portant evidences for the formation of QGP in relativistic nuclear collisions. Below 
we summarize some of the main results. 

From the comparison of relativistic hydrodynamics calculations with the 
anisotropic flow measurements, we have found that the values of shear viscosity 
to entropy ratio rj/s is on average larger at the LHC than that at RHIC; this indi¬ 
cates that the produced QGP is less strongly-coupled at the LHC. Event-by-event 
fluctuation and correlation observables associated with anisotropic flow can be uti¬ 
lized to probe the linear as well as non-linear nature of hydrodynamic response, 
and can provide tight constraints on the features of the initial conditions and final 
state correlations of heavy-ion collisions. The studies of longitudinal fluctuations, 
pre-equilibrium dynamics as well as the collective behavior in small collision sys¬ 
tems are important for understanding the anisotropic flow phenomena in heavy-ion 
collisions. Utilizing the constraint from jet quenching measurements, various phe¬ 
nomenological studies have shown that the average value of the scaled jet quenching 
parameter q/T^ is smaller at the LHC than that at RHIC; this suggests that jet- 
medium interaction is weaker at the LHC. The energy loss of full jets and the nuclear 
modification of full jet substructures can provide us more differential and detailed 
information about jet-medium interaction. The medium response to jet transport 
is important for more comprehensive understanding of jet-medium interaction and 
many jet-associated observables. 

One important future task in relativistic heavy-ion collisions is to establish a 
general theoretical framework which not only incorporates realistic hydrodynamics 
models and jet energy loss/deposition calculations, but also allows us to simulate 
the bulk matter evolution and hard jet transport together at the same time. Given 
a wealth of progresses made on the study of bulk matter and hard jets in recent 
years and many experimental data to come, it is hopeful that in the next few years 
many remaining questions can be resolved and more quantitative understanding of 
the transport properties of the QGP will be achieved. 


Acknowledgements 

This work was supported in part by the Natural Science Foundation of Ghina under 
grant No. 11375072. 



February 10, 2015 1:34 WSPC/INSTRUCTION FILE 


ws-ijmpe 


32 Guang- You Qin 

References 

1. Wuppertal-Budapest Collaboration Collaboration (S. Borsanyi et al), JHEP 1009 
(2010) 073, arXiv: 1005.3508. 

2. A. Bazavov, T. Bhattacharya, M. Cheng, C. DeTar, H. Ding et al., Phys.Rev. D85 
(2012) 054503, arXiv:1111.1710. 

3. T. Bhattacharya, M. I. Bnchoff, N. H. Christ, H.-T. Ding, R. Gupta et al, 
Phys.Rev.Lett. 113 (2014) 082001, arXiv:1402.5175. 

4. U. Heinz and R. Snellings, Ann.Rev.Nucl.Part.Sci. 63 (2013) 123, arXiv:1301.2826. 

5. C. Gale, S. Jeon and B. Schenke, Int.J.Mod.Phys. A28 (2013) 1340011, 
arXiv: 1301.5893. 

6. H. Song (2013) arXiv:1401.0079. 

7. P. Romatschke, Int.J.Mod.Phys. E19 (2010) 1, arXiv:0902.3663. 

8. P. Huovinen, Int.J.Mod.Phys. E22 (2013) 1330029, arXiv:1311.1849. 

9. M. Gynlassy and L. McLerran, Nucl.Phys. A750 (2005) 30, arXiv:nucl-th/0405013. 

10. PHENIX Collaboration Collaboration (K. Adcox et al.), Nucl.Phys. A757 (2005) 
184, arXiv:nucl-ex/0410003. 

11. BRAHMS Collaboration Collaboration (I. Arsene et al.), Nucl.Phys. A757 (2005) 
1, arXiv:nncl-ex/0410020. 

12. B. Back, M. Baker, M. Ballintijn, D. Barton, B. Becker et al, Nucl.Phys. A757 
(2005) 28, arXiv:nucl-ex/0410022. 

13. STAR Collaboration Collaboration (J. Adams et al.), Nucl.Phys. A757 (2005) 102, 
arXiv:nucl-ex/0501009. 

14. B. Muller and J. L. Nagle, Ann.Rev.Nucl.Part.Sci. 56 (2006) 93, arXivmncl- 
th/0602029. 

15. B. Muller, J. Schukraft and B. Wyslouch, Ann.Rev.Nucl.Part.Sci. 62 (2012) 361, 
arXiv: 1202.3233. 

16. A. Majumder, B. Muller and X.-N. Wang, Phys.Rev.Lett. 99 (2007) 192301, 
arXiv:hep-ph/0703082. 

17. J.-Y. Ollitrault, Phys.Rev. D46 (1992) 229. 

18. PHENIX Collaboration Collaboration (S. Adler et al.), Phys.Rev.Lett. 91 (2003) 
182301, arXiv:nucl-ex/0305013. 

19. STAR Collaboration Collaboration (J. Adams et al.), Phys.Rev.Lett. 92 (2004) 
052302, arXiv:nucl-ex/0306007. 

20. The ALICE Collaboration Collaboration (K. Aamodt et al.), Phys.Rev.Lett. 105 
(2010) 252302, arXiv: 1011.3914. 

21. ATLAS Collaboration Collaboration (G. Aad et al.), Phys.Lett. B707 (2012) 330, 
arXiv:1108.6018. 

22. GMS Collaboration Collaboration (S. Chatrchyan et al), Phys.Rev. C87 (2013) 
014902, arXiv: 1204.1409. 

23. D. H. Rischke, S. Bernard and J. A. Maruhn, Nucl. Phys. A595 (1995) 346, 
arXiv: nucl-th/9504018. 

24. D. H. Rischke, Y. Pursun and J. A. Maruhn, Nucl.Phys. A595 (1995) 383, 
arXiv: nucl-th/9504021. 

25. P. Kolb, P. Huovinen, U. W. Heinz and H. Heiselberg, Phys.Lett. B500 (2001) 232, 
arXiv:hep-ph/0012137. 

26. D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. Lett. 86 (2001) 4783, arXiv:nucl- 
th/0011058. 

27. P. Huovinen, P. Kolb, U. W. Heinz, P. Ruuskanen and S. Voloshin, Phys.Lett. B503 
(2001) 58, arXiv:hep-ph/0101136. 

28. T. Hirano and K. Tsuda, Phys.Rev. C66 (2002) 054905, arXiv:nucl-th/0205043. 



February 10, 2015 1:34 WSPC/INSTRUCTION FILE ws-ijmpe 


Anisotropic Flow and Jet Quenching in Relativistic Nuclear Collisions 33 

29. C. Nonaka and S. A. Bass, Phys. Rev. C75 (2007) 014902, arXiv:nucl-th/0607018. 

30. P. Romatschke and U. Romatschke, Phys. Rev. Lett. 99 (2007) 172301, 
arXiv:0706.1522. 

31. K. Dusling and D. Teaney, Phys. Rev. C77 (2008) 034905, arXiv:0710.5932. 

32. D. Molnar and P. Huovinen, J.Phys. G35 (2008) 104125, arXiv:0806.1367. 

33. H. Song and U. W. Heinz, Phys.Rev. C77 (2008) 064901, arXiv:0712.3715. 

34. H. Niemi, K. J. Eskola and P. V. Ruuskanen, Phys. Rev. C79 (2009) 024903, 
arXiv:0806.1116. 

35. M. Luzum and P. Romatschke, Phys. Rev. C78 (2008) 034915, arXiv:0804.4015. 

36. B. Schenke, S. Jeon and C. Gale, Phys.Rev.Lett. 106 (2011) 042301, arXiv: 1009.3244. 

37. C. Shen, Z. Qiu, H. Song, J. Bernhard, S. Bass et al. (2014) arXiv:1409.8164. 

38. W. Israel, Annals Phys. 100 (1976) 310. 

39. W. Israel and J. Stewart, Annals Phys. 118 (1979) 341. 

40. R. Baler, P. Romatschke, D. T. Son, A. O. Starinets and M. A. Stephanov, JHEP 
0804 (2008) 100, arXiv:0712.2451. 

41. P. Romatschke, Class. Quant. Grav. 27 (2010) 025006, arXiv:0906.4787. 

42. A. El, Z. Xu and C. Greiner, Phys.Rev. C81 (2010) 041901, arXiv:0907.4500. 

43. G. Denicol, H. Niemi, E. Molnar and D. Rischke, Phys.Rev. D85 (2012) 114047, 
arXiv:1202.4551. 

44. A. Jaiswal, R. S. Bhalerao and S. Pal, Phys.Rev. C87 (2013) 021901, 
arXiv: 1302.0666. 

45. D. Bazow, U. W. Heinz and M. Strickland, Phys.Rev. COO (2014) 044908, 
arXiv:1311.6720. 

46. H. Song and U. W. Heinz, Phys.Rev. C81 (2010) 024905, arXiv:0909.1549. 

47. F. Cooper and G. Frye, Phys.Rev. DIO (1974) 186. 

48. H. Song, S. Bass and U. W. Heinz, Phys.Rev. C89 (2014) 034919, arXiv:1311.0157. 

49. ALICE Collaboration Collaboration (B. Abelev et al.), Phys.Rev. C88 (2013) 
044910, arXiv: 1303.0737. 

50. P. Huovinen and P. Petreczky, Nucl.Phys. A837 (2010) 26, arXiv:0912.2541. 

51. H.-J. Drescher and Y. Nara, Phys.Rev. C75 (2007) 034905, arXiv:nucl-th/0611017. 

52. H. Song, S. A. Bass, U. Heinz, T. Hirano and C. Shen, Phys.Rev.Lett. 106 (2011) 
192301, arXiv:1011.2783. 

53. M. Miller and R. Snellings (2003) arXiv:nucl-ex/0312008. 

54. STAR Collaboration Collaboration (J. Adams et al ), Phys.Rev. C72 (2005) 014904, 
arXiv: nucl-ex/0409033. 

55. J.-Y. Ollitrault, A. M. Poskanzer and S. A. Voloshin, Phys.Rev. C80 (2009) 014904, 
arXiv:0904.2315. 

56. M. Gyulassy, D. H. Rischke and B. Zhang, Nucl.Phys. A613 (1997) 397, arXiv:nucl- 
th/9609030. 

57. G. Aguiar, Y. Hama, T. Kodama and T. Osada, Nucl.Phys. A698 (2002) 639, 
arXiv:hep-ph/0106266. 

58. R. Andrade, F. Grassi, Y. Hama, T. Kodama and J. Socolowski, O., Phys.Rev.Lett. 
97 (2006) 202302, arXiv:nucl-th/0608067. 

59. W. Broniowski, P. Bozek and M. Rybczynski, Phys.Rev. C76 (2007) 054905, 
arXiv:0706.4266. 

60. R. Andrade, F. Grassi, Y. Hama, T. Kodama and W. Qian, Phys.Rev.Lett. 101 
(2008) 112301, arXiv:0805.0018. 

61. T. Hirano and Y. Nara, Phys.Rev. C79 (2009) 064904, arXiv:0904.4080. 

62. B. Alver and G. Roland, Phys. Rev. C81 (2010) 054905, arXiv: 1003.0194. 

63. B. H. Alver, C. Gombeaud, M. Luzum and J.-Y. Ollitrault (2010) arXiv: 1007.5469. 



February 10, 2015 1:34 WSPC/INSTRUCTION FILE ws-ijmpe 


34 Guang-You Qin 

64. H. Petersen, G.-Y. Qin, S. A. Bass and B. Muller, Phys.Rev. C82 (2010) 041901, 
arXiv: 1008.0625. 

65. P. Staig and E. Shuryak, Phys.Rev. C84 (2011) 034908, arXiv:1008.3139. 

66. G.-Y. Qin, H. Petersen, S. A. Bass and B. Muller, Phys.Rev. C82 (2010) 064903, 
arXiv: 1009.1847. 

67. G.-L. Ma and X.-N. Wang, Phys.Rev.Lett. 106 (2011) 162301, arXiv:1011.5249. 

68. J. Xu and G. M. Ko, Phys.Rev. C83 (2011) 021903, arXiv: 1011.3750. 

69. D. Teaney and L. Yan, Phys.Rev. C83 (2011) 064904, arXiv:1010.1876. 

70. Z. Qiu and U. W. Heinz, Phys.Rev. C84 (2011) 024911, arXiv: 1104.0650. 

71. R. S. Bhalerao, M. Luzum and J.-Y. Ollitrault (2011) arXiv: 1104.4740. 

72. S. Floerchinger and U. A. Wiedemann (2011) arXiv:1108.5535. 

73. L. Pang, Q. Wang and X.-N. Wang, Phys.Rev. C86 (2012) 024911, arXiv:1205.5019. 

74. B. Schenke, P. Tribedy and R. Venugopalan, Phys.Rev.Lett. 108 (2012) 252301, 
arXiv: 1202.6646. 

75. C. E. Goleman-Smith and B. Mller, Phys.Rev. D89 (2014) 025019, arXiv:1307.5911. 

76. M. Luzum and H. Petersen, J.Phys. G41 (2014) 063102, arXiv:1312.5503. 

77. H. Niemi, G. Denicol, H. Holopainen and P. Huovinen, Phys.Rev. C87 (2013) 054901, 
arXiv:1212.1008. 

78. ATLAS Collaboration Collaboration (G. Aad et al.), Phys.Rev. C86 (2012) 014907, 
arXiv: 1203.3087. 

79. M. Luzum and J.-Y. Ollitrault, Nucl.Phys. A904-905 (2013) 377c, arXiv:1210.6010. 

80. F. G. Gardim, F. Grassi, Y. Hama, M. Luzum and J.-Y. Ollitrault, Phys.Rev. C83 
(2011) 064901, arXiv: 1103.4605. 

81. PHENIX Collaboration Collaboration (A. Adare et al.), Phys.Rev.Lett. 107 (2011) 
252301, arXiv: 1105.3928. 

82. STAR Collaboration Collaboration (Y. Pandit), Nucl.Phys. A904-905 (2013) 357c, 
arXiv:1210.5315. 

83. H. Niemi, G. S. Denicol, P. Huovinen, E. Molnar and D. H. Rischke, Phys.Rev.Lett. 
106 (2011) 212302, arXiv:1101.2442. 

84. ATLAS Collaboration Collaboration (G. Aad et al.), JHEP 1311 (2013) 183, 
arXiv: 1305.2942. 

85. ALICE Collaboration (A. R. Timmins), J.Phys.Conf.Ser. 446 (2013) 012031, 
arXiv: 1301.6084. 

86. C. Gale, S. Jeon, B. Schenke, P. Tribedy and R. Venugopalan, Phys.Rev.Lett. 110 
(2013) 012302, arXiv: 1209.6330. 

87. T. Renk and H. Niemi (2014) arXiv:1401.2069. 

88. Z. Qiu and U. Heinz, Phys.Lett. B717 (2012) 261, arXiv: 1208.1200. 

89. Measurement of the correlation between elliptic flow and higher-order flow harmonics 
in Pb-hPb collisions at .,fsfff=2.76 TeV, Tech. Rep. ATLAS-CONF-2014-022, CERN 
(May 2014). 

90. ATLAS Collaboration Collaboration (J. Jia), Nucl.Phys. A910-911 (2013) 276, 
arXiv: 1208.1427. 

91. R. S. Bhalerao, M. Luzum and J.-Y. Ollitrault, Phys.Rev. C84 (2011) 054901, 
arXiv:1107.5485. 

92. G.-Y. Qin and B. Muller, Phys.Rev. C85 (2012) 061901, arXiv:1109.5961. 

93. J. Jia and D. Teaney, Eur.Phys.J. C73 (2013) 2558, arXiv: 1205.3585. 

94. A. M. Poskanzer and S. Voloshin, Phys.Rev. C58 (1998) 1671, arXiv:nucl- 
ex/9805001. 

95. J.-Y. Ollitrault, Phys.Rev. D48 (1993) 1132, arXiv:hep-ph/9303247. 

96. H. Petersen, V. Bhattacharya, S. A. Bass and C. Greiner, Phys.Rev. C84 (2011) 



February 10, 2015 1:34 WSPC/INSTRUCTION FILE ws-ijmpe 


Anisotropic Flow and Jet Quenching in Relativistic Nuclear Collisions 35 
054908, arXiv: 1105.0340. 

97. K. Xiao, F. Liu and F. Wang, Phys.Rev. C87 (2013) 011901, arXiv:1208.1195. 

98. J. Jia and P. Huo (2014) arXiv:1403.6077. 

99. J. Jia and P. Huo, Phys.Rev. C90 (2014) 034905, arXiv: 1402.6680. 

100. L.-G. Pang, G.-Y. Qin, V. Roy, X.-N. Wang and G.-L. Ma (2014) arXiv: 1410.8690. 

101. ALIGE Gollaboration Gollaboration (B. Abelev et al.), Phys.Lett. B719 (2013) 29, 
arXiv:1212.2001. 

102. ATLAS Gollaboration Gollaboration (G. Aad et al.), Phys.Rev.Lett. 110 (2013) 
182302, arXiv:1212.5198. 

103. CMS Collaboration Collaboration (S. Chatrchyan et al.) (2013) arXiv: 1305.0609. 

104. PHENIX Collaboration Collaboration (A. Adare et al.), Phys.Rev.Lett. (2013) 
arXiv: 1303.1794. 

105. CMS Collaboration Collaboration (V. Khachatryan et al.), JHEP 1009 (2010) 091, 
arXiv: 1009.4122. 

106. ALICE Collaboration Collaboration (B. B. Abelev et al.), Phys.Lett. B726 (2013) 
164, arXiv: 1307.3237. 

107. P. Bozek and W. Broniowski (2013) arXiv:1304.3044. 

108. A. Bzdak, B. Schenke, P. Tribedy and R. Venugopalan (2013) arXiv: 1304.3403. 

109. G.-Y. Qin and B. Mller, Phys.Rev. C89 (2014) 044902, arXiv: 1306.3439. 

110. B. Schenke and R. Venugopalan, Phys.Rev.Lett. 113 (2014) 102301, arXiv: 1405.3605. 

111. P. Bozek, W. Broniowski and G. Torrieri, Phys.Rev.Lett. Ill (2013) 172303, 
arXiv: 1307.5060. 

112. K. Werner, M. Bleicher, B. Guiot, I. Karpenko and T. Pierog, Phys.Rev.Lett. 112 
(2014) 232301, arXiv: 1307.4379. 

113. G.-L. Ma and A. Bzdak, Phys.Lett. B739 (2014) 209, arXiv:1404.4129. 

114. Z.-W. Lin, C. M. Ko, B.-A. Li, B. Zhang and S. Pal, Phys.Rev. C72 (2005) 064901, 
arXiv:nucl-th/0411110. 

115. K. Dusling and R. Venugopalan, Phys.Rev. D87 (2013) 094034, arXiv:1302.7018. 

116. J. Nagle, A. Adare, S. Beckman, T. Koblesky, J. O. Koop et al., Phys.Rev.Lett. 113 
(2014) 112301, arXiv: 1312.4565. 

117. P. Bozek and W. Broniowski, Phys.Lett. B739 (2014) 308, arXiv:1409.2160. 

118. J. D. Bjorken FERMILAB-PUB-82-059-THY. 

119. M. Gyulassy and M. Plumer, Phys.Lett. B243 (1990) 432. 

120. X.-N. Wang and M. Gyulassy, Phys.Rev.Lett. 68 (1992) 1480. 

121. PHENIX Collaboration (K. Adcox et al), Phys. Rev. Lett. 88 (2002) 022301, 
arXiv:nucl-ex/0109003. 

122. STAR Collaboration (C. Adler et al.), Phys. Rev. Lett. 89 (2002) 202301, arXiv:nucl- 
ex/0206011. 

123. ALICE Collaboration Collaboration (K. Aamodt et al ), Phys.Lett. B696 (2011) 30, 
arXiv:1012.1004. 

124. PHENIX Collaboration Collaboration (D. Sharma), J.Phys. G38 (2011) 124082. 

125. CMS Collaboration Collaboration (S. Chatrchyan et al.), Eur.Phys.J. C72 (2012) 
1945, arXiv: 1202.2554. 

126. STAR Collaboration Collaboration (C. Adler et al.), Phys.Rev.Lett. 90 (2003) 
082302, arXiv:nucl-ex/0210033. 

127. STAR Collaboration Collaboration (J. Adams et al.), Phys.Rev.Lett. 97 (2006) 
162301, arXiv:nucl-ex/0604018. 

128. PHENIX Collaboration Collaboration (A. Adare et al.), Phys.Rev. C80 (2009) 
024908, arXiv:0903.3399. 

129. STAR Collaboration Collaboration (B. Abelev et al ), Phys.Rev. C82 (2010) 034909, 



February 10, 2015 1:34 WSPC/INSTRUCTION FILE 


ws-ijmpe 


36 Guang-You Qin 
arXiv:0912.1871. 

130. Atlas Collaboration Collaboration (G. Aad et al ), Phys.Rev.Lett. 105 (2010) 252303, 
arXiv:1011.6182. 

131. CMS Collaboration Collaboration (S. Chatrchyan et al.), Phys.Rev. C84 (2011) 
024906, arXiv: 1102.1957. 

132. CMS Collaboration Collaboration (S. Chatrchyan et al.) (2012) arXiv: 1205.0206. 

133. ATLAS Collaboration Collaboration (G. Aad et al.) (2014) arXiv: 1411.2357. 

134. CMS Collaboration Collaboration (S. Chatrchyan et al.), JHEP 1210 (2012) 087, 
arXiv: 1205.5872. 

135. CMS Collaboration Collaboration (S. Chatrchyan et al.), Phys.Lett. B730 (2014) 
243, arXiv: 1310.0878. 

136. ATLAS Collaboration Collaboration (G. Aad et al.) (2014) arXiv: 1406.2979. 

137. B. Jager, A. Schafer, M. Stratmann and W. Vogelsang, Phys.Rev. D67 (2003) 
054005, arXiv:hep-ph/0211007. 

138. Y. L. Dokshitzer, Sov.Phys. JETP 46 (1977) 641. 

139. V. Gribov and L. Lipatov, Sov. J.Nucl.Phys. 15 (1972) 438. 

140. G. Altarelli and G. Parisi, Nucl.Phys. B126 (1977) 298. 

141. PHENIX Gollaboration Collaboration (A. Adare et al.), Phys.Rev. D76 (2007) 
051106, arXiv:0704.3599. 

142. K. Eskola, H. Paukknnen and C. Salgado, JHEP 0904 (2009) 065, arXiv:0902.4154. 

143. M. Hirai, S. Kumano and T.-H. Nagai, Phys.Rev. C76 (2007) 065207, 
arXiv:0709.3038. 

144. D. de Florian and R. Sassot, Phys.Rev. D69 (2004) 074028, arXiv:hep-ph/0311227. 

145. R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl.Phys. B483 
(1997) 291, arXiv:hep-ph/9607355. 

146. R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl.Phys. B484 
(1997) 265, arXiv:hep-ph/9608322. 

147. B. Zakharov, JETP Lett. 63 (1996) 952, arXiv:hep-ph/9607440. 

148. M. Gynlassy, P. Levai and I. Vitev, Nucl.Phys. B571 (2000) 197, arXiv:hep- 
ph/9907461. 

149. M. Gyulassy, P. Levai and 1. Vitev, Phys.Rev.Lett. 85 (2000) 5535, arXiv:nncl- 
th/0005032. 

150. M. Djordjevic and U. W. Heinz, Phys.Rev.Lett. 101 (2008) 022302, arXiv:0802.1230. 

151. U. A. Wiedemann, Nucl.Phys. B588 (2000) 303, arXiv:hep-ph/0005129. 

152. U. A. Wiedemann, Nucl.Phys. A690 (2001) 731, arXiv:hep-ph/0008241. 

153. X.-f. Gno and X.-N. Wang, Phys.Rev.Lett. 85 (2000) 3591, arXiv:hep-ph/0005044. 

154. X.-N. Wang and X.-f. Gno, Nucl.Phys. A696 (2001) 788, arXiv:hep-ph/0102230. 

155. A. Majnmder, Phys.Rev. D85 (2012) 014023, arXiv:0912.2987. 

156. P. B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0111 (2001) 057, arXiv:hep- 
ph/0109064. 

157. P. B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0206 (2002) 030, arXiv:hep- 
ph/0204343. 

158. S. Caron-Huot and C. Gale, Phys.Rev. C82 (2010) 064902, arXiv: 1006.2379. 

159. N. Armesto, B. Cole, C. Gale, W. A. Horowitz, P. Jacobs et al. (2011) 
arXiv:1106.1106. 

160. L. Landan and 1. Pomeranchuk, Dokl.Akad.Nauk Ser.Fiz. 92 (1953) 535. 

161. A. B. Migdal, Phys.Rev. 103 (1956) 1811. 

162. C. A. Salgado and U. A. Wiedemann, Phys.Rev. D68 (2003) 014008, arXiv:hep- 
ph/0302184. 

163. T. Renk, J. Ruppert, C. Nonaka and S. A. Bass, Phys.Rev. C75 (2007) 031902, 



February 10, 2015 1:34 WSPC/INSTRUCTION FILE ws-ijmpe 


Anisotropic Flow and Jet Quenching in Relativistic Nuclear Collisions 37 
arXiv:nucl-th/0611027. 

164. S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, Nucl. Phys. A784 (2007) 
426, arXiv:nucl-th/0512076. 

165. S. Jeon and G. D. Moore, Phys.Rev. C71 (2005) 034901, arXiv:hep-ph/0309332. 

166. G.-Y. Qin, J. Ruppert, S. Turbide, G. Gale, C. Nonaka et al., Phys.Rev. C76 (2007) 
064907, arXiv:0705.2575. 

167. G.-Y. Qin, J. Ruppert, G. Gale, S. Jeon, G. D. Moore and M. G. Mustafa, Phys. 
Rev. Lett. 100 (2008) 072301, arXiv:0710.0605. 

168. A. Majumder (2009) arXiv:0901.4516. 

169. G.-Y. Qin and A. Majumder, Phys.Rev.Lett. 105 (2010) 262301, arXiv:0910.3016. 

170. A. Majumder and G. Shen, Phys.Rev.Lett. 109 (2012) 202301, arXiv:1103.0809. 

171. N.-B. Ghang, W.-T. Deng and X.-N. Wang, Phys.Rev. C89 (2014) 034911, 
arXiv:1401.5109. 

172. S. Gao, G.-Y. Qin and S. A. Bass, Phys.Rev. C88 (2013) 044907, arXiv:1308.0617. 

173. B. Schenke, G. Gale and G.-Y. Qin, Phys.Rev. C79 (2009) 054908, arXiv:0901.3498. 

174. G. D. Moore and D. Teaney, Phys.Rev. C71 (2005) 064904, arXiv:liep-pli/0412346. 

175. M. G. Mustafa, Phys.Rev. C72 (2005) 014905, arXiv:hep-ph/0412402. 

176. Y. L. Dokshitzer and D. Kharzeev, Phys.Lett. B519 (2001) 199, arXiv:hep- 

pli/0106202. 

177. S. Gao, G.-Y. Qin, S. A. Bass and B. Mller, Nucl.Phys. A904-905 (2013) 653c, 
arXiv: 1209.5410. 

178. J. Uphoff, O. Fochler, Z. Xu and C. Greiner, Phys.Lett. B717 (2012) 430, 
arXiv: 1205.4945. 

179. R. Abir, U. Jamil, M. G. Mustafa and D. K. Srivastava, Phys.Lett. B715 (2012) 183, 
arXiv:1203.5221. 

180. M. Nahrgang, J. Aichelin, S. Bass, P. B. Gossiaux and K. Werner (2014) 
arXiv: 1410.5396. 

181. G.-Y. Qin and B. Muller, Phys.Rev.Lett. 106 (2011) 162302, arXiv:1012.5280. 

182. G. Y. Qin, A. Majumder, H. Song and U. Heinz, Phys. Rev. Lett. 103 (2009) 152303, 
arXiv:0903.2255. 

183. R. B. Neufeld and B. Muller, Phys. Rev. Lett. 103 (2009) 042301, arXiv:0902.2950. 

184. S. A. Bass et al, Phys. Rev. C79 (2009) 024901, arXiv:0808.0908. 

185. A. Majumder and B. Muller, Phys. Rev. C77 (2008) 054903, arXiv:0705.1147. 

186. G.-Y. Qin and A. Majumder, Phys.Rev. C87 (2013) 024909, arXiv:1205.5741. 

187. JET Gollaboration (K. M. Burke et al.), Phys.Rev. COO (2014) 014909, 
arXiv:1312.5003. 

188. G.-Y. Qin, J. Ruppert, G. Gale, S. Jeon and G. D. Moore, Phys.Rev. C80 (2009) 
054909, arXiv:0906.3280. 

189. C. Young, B. Schenke, S. Jeon and G. Gale, Phys.Rev. C86 (2012) 034905, 
arXiv: 1111.0647. 

190. X.-F. Chen, T. Hirano, E. Wang, X.-N. Wang and H. Zhang, Phys.Rev. C84 (2011) 
034902, arXiv:1102.5614. 

191. J. Xu, A. Buzzatti and M. Gyulassy (2014) arXiv: 1402.2956. 

192. A. Majumder, E. Wang and X.-N. Wang, Phys.Rev.Lett. 99 (2007) 152301, 
arXiv:nucl-th/0412061. 

193. Z.-q. Zhang, D.-f. Hon and H.-c. Ren, JHEP 1301 (2013) 032, arXiv:1210.5187. 

194. T. Renk, Phys.Rev. C85 (2012) 044903, arXiv:1112.2503. 

195. G.-Y. Qin and A. Majumder (2014) arXiv: 1411.5642. 

196. C. Coleman-Smith and B. Muller (2012) arXiv: 1209.3328. 

197. B. Wu, JHEP 1110 (2011) 029, arXiv:1102.0388. 



February 10, 2015 1:34 WSPC/INSTRUCTION FILE ws-ijmpe 


38 Guang-You Qin 

198. T. Liou, A. Mueller and B. Wu, Nucl.Phys. A916 (2013) 102, arXiv: 1304.7677. 

199. E. lancu (2014) arXiv: 1403.1996. 

200. J.-P. Blaizot and Y. Mehtar-Tani (2014) arXiv: 1403.2323. 

201. Z.-B. Kang, E. Wang, X.-N. Wang and H. Xing, Phys.Rev.Lett. 112 (2014) 102001, 
arXiv:1310.6759. 

202. Z.-B. Kang, E. Wang, X.-N. Wang and H. Xing (2014) arXiv:1409.1315. 

203. A. Majumder (2012) arXiv:1202.5295. 

204. M. Panero, K. Rummukainen and A. SchOlfer, Phys.Rev.Lett. 112 (2014) 162001, 
arXiv: 1307.5850. 

205. S. Caron-Huot, Phys.Rev. D79 (2009) 065039, arXiv:0811.1603. 

206. M. Laine and A. Rothkopf (2013) arXiv:1310.2413. 

207. F. D’Eramo, M. Lekaveckas, H. Liu and K. Rajagopal, JHEP 1305 (2013) 031, 
arXiv:1211.1922. 

208. STAR Collaboration Collaboration (M. Ploskon), Nucl.Phys. A830 (2009) 255C, 
arXiv:0908.1799. 

209. PHENIX Collaboration Collaboration (Y.-S. Lai), Nucl.Phys. A830 (2009) 251C, 
arXiv:0907.4725. 

210. I. Vitev and B.-W. Zhang, Phys.Rev.Lett. 104 (2010) 132001, arXiv:0910.1090. 

211. W. Dai, 1. Vitev and B.-W. Zhang (2012) arXiv: 1207.5177. 

212. C. Young, B. Schenke, S. Jeon and C. Gale, Phys.Rev. C84 (2011) 024907, 
arXiv:1103.5769. 

213. Y. He, 1. Vitev and B.-W. Zhang, Phys.Lett. B713 (2012) 224, arXiv:1105.2566. 

214. T. Renk, Phys.Rev. C85 (2012) 064908, arXiv:1202.4579. 

215. K. C. Zapp, F. Krauss and U. A. Wiedemann, JHEP 1303 (2013) 080, 
arXiv:1212.1599. 

216. G.-Y. Qin, Eur.Phys.J. C74 (2014) 2959, arXiv:1210.6610. 

217. X.-N. Wang and Y. Zhu, Phys.Rev.Lett. Ill (2013) 062301, arXiv: 1302.5874. 

218. J. Casalderrey-Solana, J. G. Milhano and U. A. Wiedemann, J.Phys.G G38 (2011) 
035006, arXiv: 1012.0745. 

219. J. Casalderrey-Solana, Y. Mehtar-Tani, C. A. Salgado and K. Tywoniuk, Phys.Lett. 
B725 (2013) 357, arXiv:1210.7765. 

220. J.-P. Blaizot, E. lancu and Y. Mehtar-Tani, Phys.Rev.Lett. Ill (2013) 052001, 
arXiv:1301.6102. 

221. E. lancu (2013) 251, arXiv: 1305.7023. 

222. 1. Vitev, S. Wicks and B.-W. Zhang, JHEP 0811 (2008) 093, arXiv:0810.2807. 

223. A. Majumder, Phys.Rev. C88 (2013) 014909, arXiv:1301.5323. 

224. R. Prez-Ramos and T. Renk (2014) arXiv:1411.1983. 

225. G.-L. Ma, Phys.Rev. C88 (2013) 021902, arXiv: 1306.1306. 

226. J. Casalderrey-Solana, E. Shuryak and D. Teaney, J.Phys.Conf.Ser. 27 (2005) 22, 
arXiv:hep-ph/0411315. 

227. J. Ruppert and B. Muller, Phys.Lett. B618 (2005) 123, arXiv:hep-ph/0503158. 

228. T. Renk and J. Ruppert, Phys.Rev. C73 (2006) 011901, arXiv:hep-ph/0509036. 

229. I. Bouras, B. Betz, Z. Xu and C. Greiner, Phys.Rev. COO (2014) 024904, 
arXiv: 1401.3019. 

230. S. Floerchinger and K. C. Zapp (2014) arXiv:1407.1782. 

231. R. Neufeld, Phys.Rev. C79 (2009) 054909, arXiv:0807.2996. 

232. R. Neufeld and 1. Vitev, Phys.Rev. C86 (2012) 024905, arXiv: 1105.2067. 

233. T. Renk, Phys.Rev. C88 (2013) 044905, arXiv:1306.2739. 



